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Abstract. In this work we study the influence of the newtonian noise on atom 
interferometers applied to the detection of gravitational waves, and we compute 
the resulting limits to the sensitivity in two different configurations: a single atom 
interferometer, or a pair of atom interferometers operated in a differential configuration. 
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1. Introduction 

The direct detection of Gravitational Waves is one of the most exciting challenges of 
current scientific research. The first generation of ground-based optical interferonietric 
detectors, including Virgo [Tj and GEO600 [3] in Europe, and the LIGO [2] 
interferometers in USA, achieved design sensitivity and carried out several science runs, 
which set interesting upper limits on several classes of astrophysical sources |11 El El [7] . 
The construction of a "second generation" of optical interferometers. Advanced LIGO [8] 
and Virgo [9], and the new Japanese detector KAGRA [10], is well underway; thanks to 
the implementation of several technical upgrades, the advanced detectors are expected to 
come on line with a sensitivity about ten times better than first generation instruments. 
In the meanwhile, the conceptual design of third generation detectors, like the Einstein 
Telescope |1H [12], has started. For all these optical ground based detectors the 
sensitivity in the low frequency band, below 10 Hz, is ultimately limited by the so 
called "gravity gradient", or Newtonian Noise (NN) |13[ [T^. whose source is the direct 
coupling of the test masses with any mass-density change in the environment, especially 
of seismic or atmospheric origin. Atom interferometers (see [15] for a review) have been 
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proposed recently as GW detectors [IHl [UJ [TBI ESI EHl 121], on the basis of previous 
general ideas |22]. The question arises then, if the "low frequency wall" due to NN is 
relevant also for these new proposed detectors. In this paper we consider only the NN of 
seismic origin and we carry out a detailed calculation of its contribution to the sensitivity 
curve of a simple atom interferometer both in the "single detector" configuration and 
in the "coupled differential" configuration. We anticipate our conclusions: the atom 
interferometers are subject to NN in a degree similar to optical interferometers, and 
therefore will require appropriate technical solutions to overcome this limit in the 
frequency band below lOHz. The paper is organized as follows: in section [2] we consider 
a definite atom interferometer and we compute its response to a fluctuating gravity 
field; in section [3] we apply the formulas to the case of a single detector deriving the 
limits on sensitivity; finally in section H] we consider two atom interferometers operated 
in differential configurations. 

2. Newtonian noise of seismic origin in atom interferometers 

In optical interferometric GW detectors the test masses are suspended mirrors: a 
pendular suspension is indeed the best approximation on Earth for a freely falling test 
mass. In atom interferometers instead the role of test masses is played by atoms in 
free fall, hence our intent is to determine the influence of the Newtonian coupling to 
an external, time-varying mass distribution, on freely falling masses. Some general 
considerations are possible: if the effect originates from seismic noise, it is driven by 
an external masses displacement field generally of the form W (u) ~ w~^, mediated 
by a transfer function from the seism to the test masses motion behaving also as 
w^^ |T4l EH [2l], where u is the angular frequency. Therefore the effect on test masses is 
expected to be of the form 6 {u) ru~^, where 6 {u) is a kind of reduced transfer function, 
depending on the "detection device", F is a scale factor depending on the model of seismic 
waves (it is recognized that the role of main source is played by Rayleigh surface waves, 
especially the fundamental mode and few overtones |23[ [21]). To derive the actual 
expression of 6 [u) for NN in an atom interferometer, we use the ABCD formalism for 
matter waves, described elsewhere in detail [2^ [25] . 

Assume that the Hamiltonian of the motion is at most quadratic in momentum and 
position operators 
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where Pn{r) and q^ij-) are vectors of momentum and position, respectively, whereas 
a, /9, 7, 5 are suitable square matrices (note that 5 = —a^, with T indicates the 
transposed matrix), and M is the atom rest mass. Consider an atoms' beam (a 
Gaussian packet under paraxial approximation pU | [25 | [2B | [27 [ [2S ] ) which is divided and 
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recombined through a sequence of R hght-field beam sphtters, supphed by the same 
laser: from the first beam sphtter to the last one (the output port) we may identify 
two paths, conventionally labeled s and i. By exploiting the ttt theorem |25) for the 
atoms/beam splitter interactions, and the mid-point property of Gaussian beams |29]. 
the phase difference at the output port of the interferometer can be written as: 
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where kgi^i^j is the momentum transferred to the atoms by the j-th beam splitter along 
the s {i) arm, ujs{i)j is the angular frequency of the laser beam and Os(i)j is the phase of 
the laser beam at the j-th interaction, qs{i)j is the distance of j-th interaction point from 
the laser source; equal masses are assumed for the atoms along the s and i paths. The 
expression in ([2]) is manifestly gauge-invariant |20[ [25] . and the evolution of the wave 
packets can be obtained, by means of the Ehrenfest theorem, from Hamilton's equations 
for the vector x{t) [2n| l25| l28] 
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here r represents the time-ordering operator, and an appropriate perturbative expansion 
can be used to evaluate the time-ordered exponential in (j6]) pOf |25| [28] . 

As a simple reference configuration let us consider a "Ramsey-Borde" atom 
interferometer, with a Mach-Zender geometry, as outlined in figure [T] [T5| 1201 [25]. In the 
following, we will also assume that the instrument is crossed by a plane GW with "+" 
polarization and amplitude h, propagating along the x^ = z axis, perpendicular to the 
plane of the interferometer; we adopt in the following a description in Fermi coordinates, 
which represent the best approximation to the Laboratory Cartesian system |3U] . 

Assuming the same "stable" frequency for the laser beams and neglecting the steady 
proper laser phases, the phase shift formula in ([2]) becomes 
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Figure 1. A simple atom interferometer with Mach-Zehnder geometry. Continuous 
horizontal lines, and the slanted dot-dashed lines, represent atom beams. Vertical 
dashed lines represent the laser beams; the bold continuous arrows represent relevant 
momentum transferred to the atoms; g and e mark the ground and excited internal 
states of the atoms; k is the transverse momentum in h units. 



Let us assume that atoms are subjected only to a fluctuating gravitational field g(f). 
Considering dl]), (E]), (jlj) and ([7]) we have 
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and we obtain 
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We are interested in the low frequency range, where the newtonian noise is expected to 
be the limiting factor on account of its u~^ shape. We will therefore assume that the 
single atom interferometer has a linear dimension smaller than the wavelength of seismic 
surface waves, which we will assume to set also the coherence length. Introducing the 
Fourier transform g (u) of the fluctuating field we can also write 
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and we assume, in the long wavelength approximation, that g [iS) is the same at any 
point of the interferometer. Therefore the solution of the Hamilton equations ([5]) can 
be written as 



(12) 



pibjt piujto^ 






(t-to) ^iut _ 1 

e o 

lUJ w^ 



(^ 



1) 



by iterating the relation in ([5]) to the four interaction points of the interferometer 
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in figure [H setting ^3 = ^2 and defining T = t4 — t^ 
phase shift at the output port of the interferometer: 
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this is the fundamental formula to estimate the effect of the fluctuating fleld g. We recall 
that k is the unperturbed wave vector of the laser beam, corresponding to the impulse 
(in units of the reduced Planck constant h) transferred to the atom at each interaction 
point. Note also that in the limit a; — )■ the expression in (fT3l) corresponds to the well 
known static result [291 [3T] . 



3. Newtonian-Noise limit on sensitivity: the single detector case 

In the weak fleld approximation, to flrst order in the amplitude h of an impinging 
gravitational wave, the phase shift at the output of the interferometer in flgure [1] has 
been already obtained in a fully covariant way [20j . Indicating with gi the unperturbed 
distance of the flrst interaction point from the laser, and with pi the unperturbed 
momentum of the atoms, just before the flrst interaction with the laser beam, we recall 
that the Fourier transform of the phase shift, as a function of the Fourier transformed 
amplitude h of the GW, can be written as 
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in which the proper laser phases have been neglected. For a single interferometer with 
the laser source close to the device, the last term can be neglected and the more relevant 
one is the dynamic term, proportional to pi; we can also neglect the recoil term -^ . 
This expression can be directly translated into a relation between linear power spectral 
densities (LPSD), that we denote by a tilde, deflned in terms of the two-point correlation 
functions as 

{g{uj)g{uj')) = 2715 {uj-uj')g^uj) (15) 

in which the angular brackets represent the statistical average. From (fT3l) and ( TT4l) we 
obtain 
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where the overall linear size L = 2Tpi/M of the interferometer in flgure [H has been 
introduced; combining the two equations, we deduce the expression 
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for the equivalent strain hj^^ induced by the fluctuating field g {u). 

For definiteness, it is useful to discuss here the scale of the g{uj) LPSD, referring to 
typical values measured at the site of the Virgo interferometers; we recall indeed that 
we are considering the effect of an external fluctuating gravity field on freely falling 
test masses, which is the same situation experienced by the test masses of optical 
interferometers |T^ [23| |21] : even though the detailed shape of the NN affecting a 
instrument like Virgo depends on the model for the seismic sources and the superficial 
Earth layers, similar results are obtained in different cases, which can be summarized 
as follows 

V 4X(a;) ^^ 1.2 X 10'^ , H^ 

Lv w m 

where Ly = 3000m is the length of Virgo arms, X (u) is the displacement LPSD for 
a single suspended mirror, and Xgeism {^) is the measured LPSD of the ground seism [32| ; 
the factor \/4 takes into account that in Virgo the noise due to the four end-station 
mirrors adds in quadrature. Considering the relation between the mirror motion and 
its acceleration, due to the fluctuating gravitational field, g (u) = u'^X (u), we obtain 

iM.^^»M^6xlO-«^:E„„„Mx5£!; (19, 

Li Z Lj hy Lj m 

we further assume that the seismic noise measured at the Virgo site is well approximated 

by [33] 

^ ' [cj/ (27rHz)]' ' ^ ' 

assuming an atom interferometer with L ~ 200m, we obtain 

^ ^I^Hz^ (21) 

L [u;/ (27rHz)]' ^ ' 

We show in figure [2] examples of the newtonian noise of f lT7|) assuming the expression 
in (12T]) for the LPSD of the fluctuating gravitational field, and setting T = 0.4 s for the 
time of flight of the atoms between the interaction points in the interferometer. 

The zeroes represent frequencies at which the atom interferometer is insensitive 
both to the gravity gradient noise and to GW. 

4. Two detectors operated in a differential configuration 

Let us now consider the second term in (IT^ . proportional to the position qi. This term, 
already introduced in a different context [31], is a sort of "clock" term which takes into 
account the influence of the GW on the laser beam, along its path from the source 
to a well defined physical point. Its role was discussed in recent papers |35| |36| [37] 
and the most relevant new property is the introduction of qi (path of laser beam) in 
place of L (path of atom beam); so, in order to improve the sensitivity, enlarging gi 
seems in principle easier than enlarging L; however, the need of measuring the distance 
from the laser, and strong requirements about coherence and stability while maintaining 
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f[Hz] 

Figure 2. The solid curve represents ttie effect of tire gravity gradient noise on a 
single atom interferometer, with the expected w~^ behaviour and zeroes corresponding 
to frequencies at which the instrument is insensitive both to the gravity gradient 
fluctuations and to gravitational waves. For comparison, the dashed curve represents 
the model NN effect on the Virgo interferometer. 



sufficient power density for tlie laser beam, are not minor issues, wliicli prevent drawing 
too optimistic conclusions about the practicality of the configuration. Anyway the 
suggestion of adopting a two-interferometers differential configuration [21] appears very 
appealing in order to overcome the first of the mentioned problems; furthermore it may 
yield a good common-modes rejection. Under the hypothesis of a common laser source 
for two identical Mach-Zehnder atom interferometers in differential configuration, for 
which the relative distance D satisfies the condition uD/c <^ 1 (with c the speed of light 
in vacuum), from (fT^ the overall difference between the two partial phase differences 
at the output ports can be formally obtained as 

A0 (u) = 2kD sin^ {uT/2) e'^'^h (u) (22) 

where D = q[^ — q[ as anticipated. Considering also f lT3|) we obtain for the 
differential configuration 
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where the difference in the right hand side requires some discussion. In a given 
frequency band, if the two fluctuating gravity flelds ^1.2 act upon sufficiently distant 
interferometers, they will be uncorrelated, and we will obtain for the LPSD simply a 
sum in quadrature 
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displaying no conceptual difference with respect to the limits obtained for optical 
interferometers with long D arms j32]. Considering instead a low-frequency, long- 
wavelength approximation, it may be appealing the situation in which, even with two 
separated interferometers, the residual correlation leads to a partial noise cancellation 
in ([23]). 

We recall that the signals ^12 (w) are assumed to be stochastic acceleration fields in 
positions 1 and 2, projected along the direction specified by the segment D as in figure [3] 



D 



rl = r2 - D 



. AM(t) 





^ r2 



Figure 3. Geometry of the detector: atom interferometers are located at positions 1 
and 2, and a fluctuating mass element is assumed at a location r-i in a frame having 
position 2 as the origin, and a z axis parallel to D. 

We further assume to model the stochastic noise in the simplest possible way, 
namely as due to uncorrelated fluctuations in the density of the material surrounding 
the detector [T^. In other words a density fluctation AM (t) will contribute to the 
acceleration field in points 1 and 2 as 

'2 '2 

^ , , GAM(t) . G'AM(t) (^ -.\ 

Considering only the component acting along the direction separating the two points 
1 and 2, we obtain 

, , GAM(t) 
^2^ = ;^cos(e) 

9i{t) = ^^^^ ^ [r cos {6) - D] (27) 



as the contribution to the fluctuation of the acceleration field due to a single mass 
element. To obtain the total fluctuation, we need now to sum over the space. 

We first assume that the space can be considered homogeneous around the 
detectors: this is the case for instance if the instruments are placed in a deep mine. 



Newtonian Noise Limit in Atom Interferometers for Gravitational Wave Detection 9 



at a depth much larger than D. We are therefore interested in the quantity 
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which should be summed over the volume. It is convenient to evaluate the spectral 
density 
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where, following again Saulson [m], we assume the sum to be extended over volume 
elements of linear size A/2, with AM fluctuating coherently inside these regions, and 
totally uncorrelated otherwise: 

(AM (w, f) AM {uj', ^)) = 2Tt6 {oo - co') AM^ {oo, r) 5^p . (30) 

We obtain therefore 

~2 . ^ 4G2 ^, fcos(e) rcos(^)-D 

If we additionally assume that the mass fluctuations do not depend on r, we can 
further simplify, obtaining 
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where we have approximated the sum with an integral, normalizing by the volume 
element of the coherent region (A/2) . If we were to retain only the flrst term, we would 
obtain the same result as in [T3], corrected for a factor 2. The integration over the 
angular functions is directly carried out, resulting in a lengthy expression: 

/4«M ^ ''t^g'^' -g(AA) (33) 
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which, as expected, displays double poles in r = and in r = D. Both divergences 

I and at \r - D\ > |. 



are artefacts, which should be regulated introducing cutoffs r > j and at |r — D| > ^ 



However, it is now necessary to distinguish two cases 
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Short wavelength If the distance D ^ \, then the integral over r gives 
and we obtain 



HiD,\) = - + 0{— In-) (34) 



72 . ^ SQeTTG^AM^M 

Long wavelength In the long wavelength approximation the integral in (133|) can be 
carried out assuming r > ^ ^ D, obtaining 

HiD,X)='-^ + 0(^] (36) 

hence 
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it seems at first surprising that the dependence on D in this case cancels out, 
whereas one could have expected to retain a dependence, which could lead to zero noise 
in the D — )■ limit case. However, we are actually in a situation in which the instrument 
is sensitive to the gradient of the gravity acceleration (see (123|) ). and therefore, barring 
other sources of noise, the sensitivity is independent on the baseline D. 

We can now use Eq.l2 of |2l] to relate the mass fluctuations with the measured 
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where po is the density of the medium. We finally obtain 



hNN{sw) (W) - Xseism ['^) (39) 
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where we have used the relation Xu = 27tcl, with cl the speed of longitudinal seismic 
waves; typically in the ground where Virgo is installed Cl — 1000 ms~^, whereas we can 
assume po — 2.7 x lO^kgm"^ as a typical value for the continental crust. Assuming 
D ~ 200m, we obtain 

I / X 2x 10-8^ / X Hz^ ^ Cl ,„ ,,,, 

hNN(sw) (W) - 5 Xseism{UJ) X -- > — ~ 5Hz (41) 

~ , , 2x10"*^^ , , Hz tu Cl ,^^ ,,^. 

hNN(lw) {(^) - Xseism (w) X -- < — - ~ IHz (42) 
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Comparing with (1181) for the gravity gradient noise affecting the Virgo interferometer, we 
see that in the short wavelength limit, represented by (HTl) . the frequency dependence (as 
expected) is the same. Instead, in the long wavelength limit (H2l) . the NN affecting the 
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Figure 4. Comparison of models of the Newtonian Noise as seen by the Virgo 
interferometer (dashed hne) or by an hypothetical pair of atom interferometers 
operated in differential configuration (continuous line). Above a few Hz, the different 
scale of the two curves reflects the different baseline of the two instruments: 3km for 
the length of Virgo arms, 200m for the distance between the atom interferometers. 



atom interferometer has a slower growth for w — )■ 0, reflecting the presence of correlated 
noise at the two stations, that partially cancels out in (l23l) . 

The NN limit to the atom interferometer sensitivity is displayed in figure HI over 
a frequency range which runs from the long to the short wavelength regimes, assuming 
as reference the same seismic noise x seism as in ( 120|) . We display for comparison also 
the NN affecting the Virgo instrument; in the high frequency regime, the pair of atom 
interferemeters behaves like a smaller Virgo, at a scale about one order of magnitude 
larger because of the shorter baseline; in the low frequency regime is evident the milder 
growth as a; — )■ 0, insufficient though to push the noise below the limit affecting Virgo, 
unless the frequency range is extended well below 10~^Hz. 



5. Conclusions 



In this work we have evaluated the effect of fluctuations of the gravity field on the 
sensitivity of atom interferometers, thus providing an estimate of the so-called newtonian 
(or gravity gradient) noise for this kind of instruments. We have seen that a mid- 
scale atom interferometer, with a baseline L ~ 200m, is subject to a noise essentially 
equivalent to the one affecting a large scale optical interferometer, as Virgo. We have also 
found that operating two small-scale atom interferometers, placed at a similar distance 
D ~ 200 m, in differential configuration, the level of g.g. noise is generally larger, by at 
most an order of magnitude in the frequency range considered. 

We conclude that, similarly to what is foreseen for future optical interferome- 
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ters [llj, operating successfully atom interferometers in the [0.1, 10] Hz frequency win- 
dow will require to mitigate the gravity gradient noise, either by choosing very quiet, 
possibly underground sites, or by devising clever noise subtraction strategies. 

We acknowledge that this study has a limitation in the model for the gravity 
fluctuations, which is approximate; however, as it has been the case for similar studies 
carried out for optical interferometers [Mf 123] , we believe that the use of more reflned 
models will change the numerical results only by small factors, which would not alter 
our conclusions. 
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